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Abstract. It has been shown that for a certain special type of quantum graphs the 
random-matrix form factor can be recovered to at least third order in the scaled time 
t using periodic-orbit theory. Two types of contributing pairs of orbits were identified, 
those which require time-reversal symmetry and those which do not. We present a 
new technique of dealing with contribution from the former type of orbits. 

The technique allows us to derive the third order term of the expansion for 
general graphs. Although the derivation is rather technical, the advantages of the 
technique are obvious: it makes the derivation tractable, it identifies explicitly the orbit 
configurations which give the correct contribution, it is more algorithmic and more 
system-independent, making possible future applications of the technique to systems 
other than quantum graphs. 
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1. Introduction 

The Bohigas-Giannoni-Schmit conjecture asserts that the spectral correlations in 
quantum systems with chaotic classical analogue fall into several universality classes, 
depending on the symmetries of the system. In particular, the spectral correlations 
of systems with time-reversal (TR) symmetry coincide with the relevant expressions 
obtained in random matrix theory for Orthogonal Ensembles of matrices. This claim 
was supported by a multitude of numerical examples j2], but for a long time the only 
theoretical advance for individual systems (i.e. without disorder average) was Berry's 
diagonal approximation [H] . The recent work of Sieber and Richter |3] and Sieber jH] has 
renewed the hope that the universality of spectral correlations can be explained within 
periodic-orbit theory. 

One of the more convenient statistics is the Fourier transform of the spectral two- 
point correlator, the form factor, whose universal expression for Orthogonal Ensembles 
is given by the formula 

K GO e(t) = 2t-t log(l + 2r) = 2r - 2r 2 + 2r 3 + 0(r 4 ) , (1) 

when r is in the range < r < 1. For a quantum chaotic system the form factor can 
be written in terms of a double sum over periodic orbits (PO) using the Gutzwiller 
trace- formula jH]. This double sum is the usual starting point for analysis, in which 
different classes of periodic orbit pairs are identified and their contribution is evaluated 
to reproduce the small-r expansion of the prediction (JTJ). 

Berry [3] calculated the form factor, neglecting all correlations between POs other 
than exact symmetries. Within this "diagonal approximation" , he obtained the leading 
order term in the r expansion. In E] it was shown, that for uniformly hyperbolic 
and time-reversal invariant billiards on surfaces with constant negative curvature the 
second-order contribution — 2r 2 is related to correlations within pairs of orbits differing 
in the orientation of one of the two loops resulting from a self-intersection of the orbit. 
The same result, but without restriction to uniformly hyperbolic dynamics, was derived 
for a large family of quantum graphs [§]. In particular, some progress has been made 
in identifying the exact requirements on the degree of "chaoticity" of the graphs. The 
further step to third order was performed in JU], where some of the classes of orbits 
were evaluated for general graphs. However for classes of orbits which explicitly required 
time-reversal symmetry, the uniform hyperbolicity had to be assumed. At this time it 
became clear that the method used in [§| would become intractable for the third order 
in general graphs and deriving the fourth order term would be quite impossible even for 
systems with uniformly hyperbolic dynamics. 

The present manuscript presents a different method of dealing with various classes 
of orbits. The underlying idea is the repeated application of the inclusion-exclusion 
principle to obtain a decomposition of orbit pairs into sets of which only a relatively 
small proportion gives nonzero contribution. This technique has several advantages: the 
derivation becomes tractable for general graphs, the orbit classes giving the universal 
contributions are identified explicitly, the application of the technique is a relatively 
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mechanical process which decreases the chance of missing a contribution and raises 
hope for a general derivation of the expansion to all orders. It is also worth mentioning 
that the technique is relatively system-independent in that it does not use any features 
specific to quantum graphs, operating on rather abstract diagrams. 

This article is organised as follows: In Section 121 we define our model and explain 
how the form factor can be expressed as a double sum over periodic orbits. In Section El 
this sum is rewritten in terms of diagrams, representing all orbits with a given number 
and topology of self-intersections. Diagrams contributing to the first three orders 
are identified. In Section 0] we explain our method by re-deriving the second order 
contribution and then proceed to apply it to obtain the third order term. 



2. Quantum graphs and periodic-orbit theory 

We consider graphs with N vertices connected by B directed bonds. A bond leading 
from vertex m to vertex I is denoted by (m,l). Since we are considering graphs with 
time-reversal invariance, it is necessary that for any bond (m, I) there exists also the 
reversed bond (l,m). We do not rule out the possibility of loops, i.e. bonds of the form 
(m, m), which are time- reversal invariant. 

The discrete quantum dynamics on a graph is defined in terms of a B x B unitary 
time-evolution operator with matrix elements 

^m'l'M — a m'm e \ Z ) 

describing the transition amplitudes from the directed bond (m, /) to (/',m'). Here the 
Kronecker delta ensures that a transition is possible only between joined bonds and cr^, m 
denotes the vertex-scattering matrix at vertex I. The phases <p m i are random variables 
distributed uniformly in [0, 2ir] and for a fixed B they define an ensemble of matrices 
which is used for averaging below. The form factor is defined at integer times 
t = 0, 1, ... by 

K^(r) = B^(\trSr) m , (3) 

where r is the scaled time r = t/B and (■ ■ •) denotes the averaging. We are interested 
in the limit of large graphs B — > oo, keeping the scaled time r fixed 

K(t) = lim K {b \t) , (4) 

B^oo 

since this is equivalent to the semiclassical limit of chaotic systems jZl E] • It is in this 
limit that the form factor is expected to assume the corresponding universal form 

The classical analogue of the quantum graph |Hl E2] is represented by a 
Markov chain on the graph, specified by the doubly stochastic matrix M of transition 
probabilities, 

M - I I 2 - IrrW |2 (K\ 

Mm'l,lm — l°m'Z,Zml ~ \ a m'm\ ■ \°) 

Matrix elements of powers of M give the classical probability to get from bond (m, /) 
to bond (k, n) in t steps 

M 



P {t) 

r {m,l)^(k,n) 



(6) 

nklm 
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Under very general assumptions it can be shown that the dynamics generated by M is 
ergodic and mixing ^3], i- e. for fixed B and t — > oo all transition probabilities become 
equal 

P (V>(M) ^ aS * "> °° V ( m ' Z )' n ) • 

However, since in (J3J) the limits B — » oo and £ — * oo are connected by fixing r, we need 
a stronger condition such as 

P^ B iUk, n) ^B- 1 as B - oo V(m, /), (*, n) . (8) 

The requirements on the speed of convergence depend on the order to which agreement 
with (JTJ) is required. To avoid unpleasant estimates we will restrict ourselves to graphs 
for which the convergence is faster than any power of B. 

Another, somewhat related, condition is that following a self-retracing orbits of 
length 0(B) has negligible probability. More rigorously, define matrix R by squaring 
the elements of the matrix M, 

d _ I q( B ) |4 _ I (0 |4 /q\ 

£*-m'l,lm — \° m 'l,lm\ ~ \ u m'm\ • \ J J 

We require that for any 7 

Y \W B ] ->Q as B ^00. (10) 

n,k,l,m 

faster than any power of B. 

The above conditions are in agreement with what is known about various special 
classes of graphs. The conditions are satisfied, for example, by complete graphs with 
either Neumann or Fourier vertices fOj where numerical evidence strongly supports 
the BGS conjecture [HI E] and are not satisfied by graphs with Dirichlet vertices and 
Neumann star graphs [E], which are known to produce non-RMT statistics. 

A connection between the quantum form factor (J3J) and the classical dynamics given 
by (jHJ) can be established by representing the form factor as a sum over (classical) POs. 
A PO is a sequence of vertices P = [p\, . . . ,p t ] defined up to a cyclic shift and such that 
(PiiPi+i) i s a bond of the graph for every i = 1, . . . , t. 

By expanding the matrix powers of S in (JHJ) and averaging over phases <f> one arrives 
at the PO expansion of the form factor El] 

k^(t)=^j:a p a* q , (ii) 

P,Q 

where Ap is the product of cr-matrices along the orbit P: Ap = cr^ 3 a^ 3 J 4 ■ ■ ■ cr^X, the 
star dentoes the complex conjugation and, most importantly, the sum is taken over only 
those pairs of orbits that visit the same set of bonds, or their time-reverses, the same 
number of times. 
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Figure 1. A schematic representation of an orbit with a self- intersection at a vertex 
a and its partner orbit (a). There are orbits, however, for which the position of the 
intersection point a is ambiguous. The pair of orbits shown on (b) can be fitted into 
the pattern (a) with either c or d playing the role of the intersection vertex. 



3. The expansion in self-intersections of the periodic orbits 

3.1. From orbits to diagrams 

The calculation of the form factor is now reduced to a combinatorial problem: ensure 
that P and Q pass through the same non-directed bonds. This can be done by composing 
P and Q from the same segments, or arcs, which would appear in P and Q in different 
order and/or orientation. We classify the pairs of orbits in the following manner. We fix 
a transformation: a permutation of arcs followed by the time reversal of selected arcs, 
and then sum over all orbits P, Q related by this transformation. The clearest way to 
represent a transformation is graphical, hence we refer to them as diagrams. The sum 
over all diagrams finally gives the form factor. 

The main problem with this approach is to ensure that each orbit pair P, Q 
is counted once and only once. This is difficult because for a given pair P, Q the 
identification of the arcs and their permutation, transforming P into Q, is not necessarily 
unique. As a simple example consider the diagrams on Fig. ^ Part (a) gives a schematic 
representation of the pair P, Q where Q is obtained from P by reversal of arc 2. Such 
orbits were considered in [Hj and were shown to give the contribution — 2r 2 . One of the 
difficult points of the derivation was correct counting of the orbits which do not merely 
cross but follow themselves for at least one bond. An example of such orbit is given in 
part (b), where the orbit crosses itself along the bond (c,d). For such orbit, there are 
two possibilities to identify arc 1: either as ^ a ^ b ^ or as ^ c ^ a ^ b ^ c 
thus taking either c or d as the intersection point. 

To avoid this double-counting, jU] imposed a restriction on arc 1. Denoting the 
first vertex of the arc by s\ and the last one by /i, it was demanded that f\. This 




Figure 2. Topology of orbits contributing to form factor at first and second orders. 
Evaluating the NTR1 orbits we recover the diagonal approximation to the form factor. 



ensured the unique choice of the arcs and the intersection point. In the example above, 
the valid choice is — > a — > b — >. 

Unfortunately, even with this restriction, there were orbits that should not be 
counted: those exceptions had self- retracing arc 2. The self-retracing arc did not change 
under time-reversal and as a result, the orbit Q was identical to the orbit P, forming 
a pair which was already counted in the diagonal approximation. The contribution of 
such orbits had to be subtracted explicitly. 

In this paper we present a different counting technique, which is easily extendable 
to more complicated diagrams, avoids the introduction of exceptions and explicitly 
identifies the orbits which give the generic contribution. We will first illustrate our 
approach by re-deriving the — 2r 2 contribution and then proceed to calculate the third 
order correction. But before we can do it, we need to introduce the notation and discuss 
some preliminary matters. 

3.2. Notation 

If we consider P as a single arc with no intersections, Q = P and Q = P are the 
only options. Here the bar denotes the operation of time-reversal. The corresponding 
diagram have a simple circular shape (the first diagram of Fig.EJ. Summation over these 
orbit pairs is nothing other than the diagonal approximation. It produces K\ = 2r. In jH] 
orbit P was treated as two arcs, 1 and 2, joined at a single intersection a, corresponding 
to an 8-shaped diagram, Fig. 121 The contributions of such orbits were found to give 
rise to the second-order term in ((TJ K 2 = —It 2 . In this paper we calculate the r 3 - 
contribution by assuming that P contains three or four arcs connected at intersections. 

We denote arcs by numbers 1,2, .. . and the intersection points by Greek letters 
a,P, . . .. An arc can be identified by a sequence of vertices, which does not include 
the intersection vertices, or, alternatively, by a sequence of bonds, which includes the 
bonds from and to the intersection points. The length of the zth arc is denoted by £j 
and is defined as the number of bonds in the arc (which is one more than the number 
of vertices in the arc). The sum of the lengths of all arcs gives t, the length of the orbit. 
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The length of an arc is at least one. Given an arc i leading from a to (3 we denote 
the first vertex following a by Sj and the last vertex before (3 by f]. In the degenerate 
case when the arc going from a to f3 is the single bond (a, (3) and does not contain any 
vertices (£j = 1) our definition implies = (3 and fi = a. 

The diagrams contributing at the third order are shown on Figs. 01 For a discussion 
why only these particular diagrams contribute at third order in r we refer the reader 
to PU]- in any diagram, the arcs forming orbit P and its partner Q are identical, but 
the way they are connected at the intersections differs. The orbit P is given by the 
connections drawn as continuous lines, while its partner orbit Q is given by connections 
drawn as dotted lines. The orbits P and Q are also written as a symbolic code to the 
right of each diagram: a path that goes from the beginning of arc 1 to vertex a then 
on arc 2 to vertex (3 and so on is denoted as la2/3- • •. The diagrams divide into two 
classes, NTR and TR. In the NTR-diagrams all the arcs of Q have the same orientation 
as the corresponding arcs in P, while in the TR-diagrams some of the arcs of Q are 
time-reversed. For a system with no time-reversal symmetry, only the NTR-diagrams 
are possible. In our case, diagrams in both classes contribute and r 3 -contributions to 
the form factor is a sum of five terms 

K% = 2 (i^NTR3a + -^NTR3b + -^TR3a + -^TR3b + -^TR3c) • (12) 

The factor of two is due to the fact that for every diagram in Fig. El there is another one 
with Q replaced by its complete time- reversal, Q, which gives an identical contribution. 

The NTR-diagrams were treated in the general case in ^Uj and it was shown that 
-^NTR3a + -^NTR3b = 0. In this manuscript we derive the contribution of TR-diagrams 
using a slightly different method. 

3.3. Avoiding double- counting 

There are degeneracies in the diagrams which can be accounted for by simple prefactors 
multiplying the contributions such as the factor of two in (fT2"|) . These factors were 
derived in ^Qj and are listed in Fig. El next to the diagrams. 

Another potential source for double-counting of orbits are tangential intersections, 
already mentioned in Section 13.11 Double-counting can be avoided using a method 
outlined in jH] and followed in [HI]: the intersection point is uniquely defined by ruling 
that if there is an ambiguity then the intersection is as far to one side as possible. As an 
example we refer to the ambiguous intersections in Fig. UJLb) and insist the intersection 
is as far to the left as possible. That is, we demand that the vertices a and b are distinct. 

Contrary to what was done in jTUj, we will not fix the restrictions for the TR 
diagrams. The restrictions we choose will depend on the lengths of individual arcs, as 
will be shown in Section |3J 

It should be emphasised that the orbits with tangential self-intersections are 
responsible for non-zero contributions to the form factor and their correct treatment 
is absolutely crucial to the derivation. 
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TR3a 




TR3b 



orbits: la 2(3 3a4|3& Ia4(3 3a2(3 
factors: 1/4 



orbits: la 2a 3 a& la 3 a 2 a 
factors: 1/3 



orbits:la2(3 3a4(3& la 4(3 2 a3 (3 
factors: 1/2 




orbits: la 2(3 3 (3 4a & la 4(3 3 (3 2a 
factors: 1/2 



TR3c 




orbits: la 2a3 a & la 3 a2 a 
factors: 1 



Figure 3. Topology of NTR3a. NTR3b. TR3a, TR3b and TR3c. In each case a pair 
of orbits is shown, one follows the solid line throughout, the second follows the solid 
line except at the intersections (denoted by circles) where it follows the dotted line. 
Each circle represents a single vertex where a self-intersection of the orbit occurs. Next 
to each topology we give the symbolic code for the pair and the corresponding weight 
factor (Section HO)! . 
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3.4- Orbit amplitudes 

Before we can attempt the summation over all orbit pairs P, Q within a given diagram, 
we need to understand the structure of the product ApA* Q appearing in (JTTJ). We 
consider the diagram TR2 as an example. Let arc 1 be of length t\, consisting of the 
vertices [s\, x 2 , £3, . . . 2^-2, fi]- Then both Ap and Aq, will contain factors cr^l, o'^ 3 2 J 1 , 
a x X 4x\ ■ ■ ■ a a%ii-l • Thus when we evaluate the product A p Aq, the contribution of the 
arc 1 will come in the form 

which is the classical probability of following arc 1 from bond (a, Si) to bond (f\,a)$. 
Analogous construction leads to the probability P 2 of following the arc 2, here we need 
to remember that the matrices a are symmetric. The factors not yet accounted for in 
Pi and P2 are the transition amplitudes picked up at the intersection vertex a: 

A P A* Q = Pl x P 2 x aflafX x (14) 

To evaluate the contribution of the diagram, (JT4*j) must be summed over all free 
parameters, namely all intersection points and all possible arcs connecting these points. 
The latter summation includes a sum over the lengths of these arcs with the restriction 
that the total length of the orbit is t. 

The summation over all the intermediate vertices £25X3, • • • > x ti-2 along arc 1 can 
be performed immediately, since it is unaffected by the restrictions discussed in the 
previous subsection. This summation adds the classical probabilities of all possible 
paths leading from bond (a, s\) to bond a) in t\ — 1 steps and results consequently 
in the classical transition probability P^~^^^ ^ given by ©• Analogous summation 
over the other arc produces P^7]\./ f „\- 

The remaining summation is over the lengths ti of all arcs, the first and the last 
vertex Sj and fi of all arcs % with ti > 1 and the intersection points like a. Here we 
try to use the fact that for sufficiently long arcs the transition probabilities can be 
replaced by B~ l according to (jHJ). Then the sum over vertices decouples into a product 
of sums associated with the intersection vertices, which can finally be evaluated using 
the unitarity of the vertex-scattering matrices a. 

4. Summation of TR diagrams 

4.1. TR2 

As was mentioned earlier, we develop a slightly different technique to deal with the TR 
orbits. We will illustrate it by first considering the TR2 contribution. In |Hj we dealt 
with the intersection point ambiguity by imposing a restriction s 2 7^ f 2 - The summation 

\ Pi = 1 if arc 1 contains no vertices, i. e. if t\ = 1. 
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then would take the form 

2B 



iW*0 = ^ (15) 



{U} 

X £ £ ^TR2 X PtR2 X AtR2 , 

Q {Si,fi} 



where 



TR2 = ^hs^hs^s 2Sx ^hh ( 16 ) 
^TR2 = ^,«i)-f(/ 1 ,a)- F (a,a 3 )-(A.a) ( 17 ) 

A TR2 = (l-5 S2/2 ). (18) 

Since at least one of the arcs must be long, we can approximate the corresponding 
a ) by its ergodic limit 1/B. If arc 1 is long, we can perform the summation 
over si (and f\) 

(i - <WJ E°jU<Ci = (i - <W 2 )<W 2 = o (19) 

Si 

to show that the contribution of such orbits is zero. However, the case when arc 2 is 
long is not nearly as easy. We can still sum over s 2 or f'2, but the restriction, which 
helped us in the first part, now stands in the way (see jU] for details). 

The idea that we are going to use is as follows: we can change the restrictions 
depending on the arc lengths. In TR2 case, if arc 1 is long, we stick with the restriction 
s 2 ^ /2, and if arc 2 is long, we switch to the restriction si ^ f\. In the first case the 
result of the summation is 0, as was shown above. In the second case we approximate 
P(a~s2)^(f2 a) by an d perform the summation over s 2 to obtain 



<WJ £ = C 1 " S s lh )S slfl = 0. (20) 



) r I J 

si 



This is not what we should be getting: the overall result should be — 2r 2 , not 0. The 
reason for getting is that there are orbits which were counted in both sums. To explain 
this, we need to define more carefully what is meant by "arc % is long" . 

We split the set of all TR2 orbit pairs into two sets. The set A will contain orbits 
satisfying t\ > t/2 and s 2 7^ / 2 - The set B will contain orbits satisfying t\ < t/2 and 
Si y£ f\. The sum of contributions coming from each set is zero, as was shown in the 
previous paragraph. However, these sets are not disjoint. Fig. 0] illustrates both sets, 
denoted by A and B, and the shape of the orbits belonging to their intersection A H B. 
These are orbits with a tangential intersection which is roughly in the middle of them: 
t\ < t/2 but t\ + 2a > t/2, where a is the length of the intersection and ti is now counted 
from the end of the tangential intersection (i.e. from vertex a ). We need to subtract 
the contribution of such orbits to recover the full term i^TR2(T)- 

To calculate the contribution of the orbits in A n B, we notice (see |Appcndix A ) 



that we get the same contribution if we modify the set A D B by setting a = 1 and 
dropping the constraints. Thus we can write for this contribution 

t 2 t +1 t +1 

ti,t 2 010,011 
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2 A 

ti > t/2 




2 B 

h < t/2 




a 



a 




AnB 
h < t/2 

h + 2a> t/2 



Figure 4. The TR2 diagram and its partition into two sets of orbits. The total 
contribution of the orbits in each set is easily shown to be zero. The result comes from 
the orbits lying in the intersection of the two sets, AnB. 



where t = t\ + t% + 2 and, because of the tangential intersection, the unitary factors 
at the intersection point got squared and absorbed into the classical probability. In 
fact, the whole expression is purely classical now. The length of the arc 1 must satisfy 
t\ < t/2 and t\ + 2 > t/2, thus there are only two possible values that ti can take, 
t/2 — 1 and t/2 — 2. The length of the arc 2 satisfies t2 =t — 2 — ti, therefore both arcs 
are long and the corresponding probabilities can be approximated by 1/B, 



\Ar\B\ 



t l 



2B 3 



E E i 



t 2 t 2 
-—2 y i = — 

2B 3 ^ B 2 



(22) 



ti,t 2 ao,ai 

where the last sum is taken over all connected pairs of vertices «o and ct\ and the factor 
of 2 appeared because there are only two possible choices of the lengths {ti, £2}- To get 
the final result for K^ 2 {'t)i we need to multiply (J22|) by two once more, to account for 
the time-reversal symmetry, and subtract it from the previous result (zero) since these 
orbits were double-counted when we calculated the contributions of the sets A and B. 
Hence we obtain the sought-after result i^TR2( r ) = — 2r 2 . It is important to note that 
in the above derivation we could have used any 0(t) value instead of t/2. 

Now our strategy of dealing with the diagrams of higher order is clear. We will 
strive to partition the orbits belonging to each diagram into several sets, such that the 
contribution of each set is easily seen to be zero. The intersections of these set will 
provide the true contributions to the form factor. Hopefully, the arcs of the orbits in 
the intersections will be long and thus the ergodic approximation can be employed to 
simplify the calculations. 



4.2. Partition of TR3a 



Before we start making partitions of the diagrams, we stress that there is no unique 
"right" way to partition. The partition that we are presenting here is only one of those 
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h + t 2 > t/2 
h < t/4 
sum over f 2 

D 

h + t 2 < t/2 
h < t/4 

sum over 



b n (C u D) 

h + t 2 + 2b> t/2 
h + t 2 < t/2 
h + b< t/4 
sum over s 2 

Figure 5. The TR3a diagram and its partition into four sets of orbits. The orbits of 
the sets A and B have identical structure but different lengths. The same applies to 
the sets C and D. A consequence of this is that Ap\B = CP\D = (i) which implies 
that all 3-way intersections are empty. The contribution of the orbits in each set can 
be easily shown to be zero. Also drawn on the Figure are three of the four 2-way 
intersections among the partition sets of the TR3a diagram. The contribution of each 
is zero. 

that we considered. All partitions produce the same results, as they should, but the 
present one seemed to be the least difficult. 

The suggested partition for the diagram TR3a is given by Fig. El (upper part). Each 
of the 4 subsets of TR3a is equipped with its own set of constraints, which are shown by 
the arrows. For example, the set A has constraints f\ ^ fa and S3 7^ f'4. The constraints 
are chosen in such a way that the contributions of each set evaluate to zero: for example, 
the arc 1 is long in the set A, so we can approximate the contribution of the arc by 
it's ergodic limit and subsequently sum over s\ to obtain zero. The constraints also 
take into account the symmetry of the diagram: under the permutation of arcs 1^3 
and 2 «-> 4, the set A is interchanged with C and B is interchanged with D. Also, it 
is immediately obvious that the sets A and B are disjoint; so are the sets C and D. 
Therefore we do not have to worry about 3-way intersections of the partition sets. 

Three of the 2-way intersections are straightforward to evaluate, see Fig. El The 





An D 

ti + t 2 + 2b> t/2 
h + t 2 < t/2 
h + b> t/4 
t 3 + b< t/4 
sum over s 4 




Form factor of quantum graphs 



13 




E = Af]C 

h + t 2 + 2b> t/2 
h + t 2 < t/2 
t\ + b> t/A 
t 3 + b> t/A 



Splits into the following sets 




+ 2x J 



t 1 + t 2 + 2b> t/2 
h + t 2 < t/2 
tx + b> t/A 
t 3 + b> t/A 

sum over sx 

E 2 

t l + t 2 + a + 2b> t/2 
h + t 2 + a< t/2 
t 2 = 
t 3 + b> t/A 



or 



E 3 

t! + t 2 + a + 2b>t/2 
h + t 2 + a< t/2 
h + b< t/A 
tx + a + b>t/A or 
t 3 + b> t/A 



U = 

tx + b> t/A 



t 3 + b< t/A 
t 3 + a + b>t/A 
tx + b> t/A 



Figure 6. The last of the 2-way intersections among the partition sets of the TR3a 
diagram. To evaluate it we represent it as the set Ei plus the orbits in the set E 2 
minus the orbits of the set E 3 . The contribution of the set E\ is zero. 



intersection AC\D has t\+t 2 m t/2, therefore t 3 + t± ~ t/2 but t 3 + b < t/A, therefore the 
arc 4 is long and we can sum over obtaining zero. The intersections Bf]C and BC\D 
can be considered together, as B H (C U D). For this set we again have t± + 1 2 ~ t/2 
and t± is small, thus enabling us to sum over s 2 . 

The last intersection, A n C = E is the hardest, see Fig. El Arcs 2 and 4 are 
potentially short and summation over any other arc is obstructed by the constraints. 
To overcome this, we change the constraint fx ^ f 3 to s 2 ^ s 4 , thus producing the set 
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Ei. There are certain orbits, however, that belong to E and do not belong to E\ and 
vice versa. These are put into the sets E 2 and E 3 , correspondingly. The orbits from the 
set Ei can be easily shown to produce zero contribution. The contribution of the sets 
E 2 and E 3 are not zero and will be evaluated separately. 

4.3. Partition of TRSb 

The diagram TR3b is split into three sets, Fig. Here r is a large number, of order 
t/4. The reason we do not just put it equal to t/4 is to highlight the interaction of the 
diagram TR3b with TR3c. 

To show that contribution of the sets A and B is zero, simply sum over Si and S3, 
correspondingly. In set C, at least one of the arcs 2 and 4 is long, thus enabling us to 
sum over s 2 or S4. It is easy to see that the sets A and B are disjoint. 

The intersection A n C is also shown on Fig. One of the arcs 2 and 4 must be 
long. To be able to sum, we need to drop the restriction f 2 7^ S4, producing the set 
Fi. By summing either over f 2 or S4 (depending on which arc is long), we show that 
the contribution of Fx is zero. However, the set F\ is larger than A H C: it includes the 
orbits sketched in the sets F 2 , those can be obtained if, say, arc 4 is one bond long and 
f 2 = S4. The set F 2 provides a non-zero contribution and will be evaluated separately. 

In a similar fashion we treat the intersection B n C, which has to be split into the 
sets F 3 , F 4 and F 5 . The contribution of the set F 3 is zero and the sets F 4 and F 5 will 
be treated separately. 

44. Partition of TR3c 

The diagram TR3c is very special as the orbits belonging to it can also be obtained 
as orbits from TR3a or/and TR3b when one of the t{ is equal to zero. Our job is to 
synchronise the partition of TR3c with the partitions we chose for TR3a and TR3b. 
The relations between TR3c and other diagrams are summarised in Table 

An immediate consequence is that all TR3c orbits with fi = f 2 were already counted 
in TR3a: no set in TR3a has the restriction s 2 7^ S4. Also some of the orbits from TR3c 
with s 2 = f 3 are counted in the sets A and C of TR3a. To understand precisely which 
orbits were counted, we notice that a TR3c orbit with s 2 = f 3 corresponds to a TR3a 
orbit with Si = f 2 and t 2 = 1. Such orbit can only belong to the set A, requiring 
t[ + t' 2 > t/2 and t[ > t/4, which is equivalent to t ± > t/2 - 2 (see the "Lengths" 
column of the table). Same applies to the reversed orbits from TR3c with s 2 = f 3 which 
were counted in TR3a set C. Thus orbits from TR3c with ti > t/2 — 2 must have the 
restriction s 2 7^ f 3 (set A of Fig. |HJ). We also note that orbits from TR3c with s\ = f 2 
and £3 > t/2 — 2 were already counted in TR3a sets B and D, warranting the use of 
corresponding restriction for the set B of Fig. |H| 

The relation between TR3c and TR3b is simpler: if the lengths of an orbit from 
TR3c satisfy ti < r and £3 < r, then we must impose restrictions Si 7^ f\ and S3 7^ f 3 , 
which is reflected in the set E of Fig. |H1 
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A 

h>r 

sum over s\ 

c 

t\ <r 
t 3 <r 

sum over s 2 or S4 



B 

h <r 

sum over S3 






Anc 



Decompose into the following sets 
3 F x 




h < r 

h + 2a>r 

t 3 <r 



h <r 

ti + 2a > r 

t 3 <r 

sum over /2 or S4 

ti < r 

ti + 2a > r 

*3 < r 

t 4 = or t 2 = 
6 > 1 




BflC ti + 2a < r 

*3 < r 

t 3 + 2b>r 



Decompose into the following sets 




h < r 

t 3 <r 

t 3 + 2b>r 

sum over jfc or S4 

ii < r 

ti + 2a > r 

*3 < r 

*3 + 26 > r 

ti < r 

*3 < r 

t 3 + 2b>r 

t 4 = or t 2 = 

a > 1 



Figure 7. The TR3b diagram and its partition into three sets of orbits. The 
contribution of the orbits in each set can be easily shown to be zero. Below the 
line we sketch the intersection of the sets A and C, re-partitioned as F± minus F2 and 
the sets B and C, re-partitioned as F3 minus F4 and F§. 
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t'i = 


tx-2 


t'2 
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t'3 = 


t 3 
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= 1 
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t'i 


= t3 


t'2 


= 1 




t'3 = 


h 


- 2 


^1 
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+ 1 


«3 = /3 




S4 = f2 




t'i 
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t'2 
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+ 1 


t'3 = 


ts 
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t'i 


= 1 








S2 = fl 




t'i = 


h-2 


t'2 


= 1 




t' 3 = 


h 




t'i 
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Table 1. Summary of conversions between degenerate orbits of TR3c and two other 
diagrams. Here i,; denotes the length of ith arc in the diagram TR3c and t\ denotes 
the same but for the diagrams TR3a and TR3b; an 'R' denotes that the diagram is 
reversed. For example, the TR3c orbit with fi = / 2 can be also considered as a TR3a 
orbit with S2 = S4 and lengths of the arcs given by t'-y = t\ — 1, t' 2 = 1, t' 3 = t% — 1 and 

= h + 1 

Bringing the above information together we fix the partition of the diagram TR3c, 
Fig. El All five partition sets are disjoint and the contributions of the sets A, B, C and 
E evaluate to zero by summing with respect to si, S3 and S3 and S2 correspondingly. 
The set A is slightly special since when £3 > t/2 — 2, there is an additional restriction 
s i 7^ $2-, but in this case one can sum over S3 to obtain the zero result. 

Despite the length restriction on the sets C and D, there are still orbits belonging to 
these sets and the set C of TR3b. These orbits are sketched on Fig. El When considering 
the intersection of C and TR3b, we assume a > since when a = 0, we can sum over 
S2 to get zero. The contributions of this set will be evaluated separately. 

The intersection of D and TR3b can be re-partitioned as Gi plus G2 minus G3 with 
the set G\ giving zero contribution (sum over f'2). We will show that the contributions 
of the sets G 2 and G3 cancel each other exactly. 

The last set to be decomposed is the set D of TR3c, see Fig. Two dashed 
restrictions are denoting the situation when S3 = f'3 and S2 = ^3 cannot be satisfied 
simultaneously . Out of five partition sets, the sets Di and D 5 are immediately evaluated 
to zero, the contributions of the sets D2 and D3 clearly cancel each other, thus only the 
contribution of the set Di survives. 
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Figure 8. The TR3c diagram and its partition into five sets of orbits. In the upper 
right corner we illustrate the transitions through the intersection vertex undertaken 
by the original (solid lines) and the partner (dashed lines) orbits. The contribution of 
all sets but D can be easily shown to be zero. The set D is re-partitioned separately 
in Fig. CHI 



4-5. Evaluating significant diagrams 

If we denote the contribution of a set A by \A\, the total contribution to the r 3 term of 
the form factor expansion takes the form 

2(#TR3a(r) + K TR3h (r) + # T R3c(t)) (23) 

= 4 (|*2| + \F*\ + \F 5 \ - \E 2 \ + \E 3 \ - 2\H\ - 2\D 4 \ - 2\G 2 \ + 2|G 3 |) • 

The sign of each contribution is the sign displayed next to it on the corresponding 
diagram inverted if this diagram is from a 2-way intersection of some sets (this is true 
for all sets but D4). The overall factor of two is due to the time-reversal symmetry and, 
in case of the sets F and E, it is cancelled by the diagram symmetry factor of 1/2, see 
Fig. El 

When evaluating the contributions we use the following rules of thumb: (1) the 
length of the tangential intersections, denoted in the diagrams by letters a and/or b, 
can be set to 1; (2) the restrictions can be ignored (see Appendix A| for a discussion of 
such approximations). We also notice that all orbits belonging to the above diagrams 
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H = C(TR3c) n C(TR3b) 

t\ < r 
t 3 <r 
t 3 + 2b>r 

(a > otherwise sum over S2) 



G = £>(TR3c) n C(TR3b) 

t\ < r 
ti + 2a>r 
h <r 

Decompose into the following sets 



h < r 

t 1 +2a>r 

t 3 <r 

sum over fi 



t\<T 
ti+2a>r 
h<r 
b > 



h <r 

ti + 2a>r 

h<r 



Figure 9. The set of orbits belonging to both set D of TR3c and set C of TR3b is 
re-partitioned as the set G\ plus G2 minus G3. The contribution of the set G\ is zero. 

give positive contribution (i.e. no more unitary factors left), which shows we are on the 
right track. 

We begin by showing that the contributions of the sets G2 and G 3 cancel each 
other. On Fig. ^2 we introduced some extra notation necessary for the calculation. The 
contribution of the set G2 can be written as 

\C I = 2 P tl+1 0*2+1 pt-t+i pi (2A) 

I 21 " / j / j 1 (a,ai)—>(ai ,o) (01,0)^(0,02) (0,02)^(02,0) (02 ,0)— >(q,oi ) ' V^V 
Ui<r cti 

where the summation is taken over as such that there exist bonds (a, ax) and (a, a?) 
and over the ranges of t^s outlined on the diagram. The two possible choices of ti are 
expressed as the factor of 2 and thus, due to the sum t\ + ti + t 3 + 4 being fixed, the 
sum is essentially over t 3 < r. The arcs 1 and 2 are (relatively) long, therefore the 
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L>(TR3c) t 3 < r 

r < t x < t/2 - 2 

Decompose into the following sets 



h 



f{— 



S3 

— h 

' S 2 



Do 



t 3 <r 
r < tx < t/2 -2 
sum over si and f\ 



t 3 <r 

t t <t/2-2 
ti + 2a> t/2 



t 3 <r 
t\<r 
ti+2a>r 



t 3 + 2a>r 

h + 2a <t/2-2 

h>r 




D 5 b>0 

t 3 + b <r 

r < t x < t/2 - 2 

sum over S2 



Figure 10. Re-partitioned set D of the TR3c diagram. Two dashed restrictions are 
denoting the situation when S3 = and S2 = cannot be satisfied simultaneously. 
The sets D\ and D5 give zero contribution and the sets D2 and D3 produce 
contributions that cancel each other. 



corresponding transition probabilities can be approximated by their ergodic limit of 
1/B. We further sum over ot\ to obtain 

\r I = — V V P t3+1 P 1 = — V V P tA+l 

I 21 D2 l—< Z^ (a,a2)^(«2,«) ("2, a)— »(a, ai) r?2 («,a2)^(o2 ,a) ^ ' 

13 tz<r oii D ts<ra,a2 



We evaluate the contribution of G 3 is a similar way, 

I 31 Z_/ / ' (o?2,ai)^(ai,02) (ai,Qf2)^(ai,02) (ai,a2)— »(o2,oi) 



*3<r ai,ct2 
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G, 



h < r 
ti + 2 > r 
h<r 



G 3 (+) ti < r 

ti + 2 > r 
t 3 < r 




# (-) ti < r 
t 3 < r 
*3 + 2 > r 




^5 (+) 



t\ < r 
t 3 <r 
h + 2>r 

u = o 



or t 2 = 




1 w P 

Figure 11. First three of the sets giving non-zero contributions 



= — V V P t3+1 (9G\ 

D2 Z_< (ai,a 2 )^(o 2 ,ai) ^ ' 

<3<rai,a2 

and notice that the final expressions for IG2I and \Gs\ are identical and thus cancel each 
other. 

The sets H and F$, Fig. ^2 on the other hand, are identical themselves. The factor 
of 2 before the contribution of H is compensated by the factor of 2 in the diagram F 5 
itself. Thus the total contribution [^5 1 — 2\H\ = 0. 

The contributions of the sets E 2 , E 3 and F 2 are best grouped together. We start 
with E 3 (evaluating only one of the two cases and multiplying it by 2), Fig. To 
simplify the notation when dealing with numbers of order t, such as t/2 and r, we 
ignore the constant corrections as they do not influence the limiting behaviour. The 
constraints imply t\ = r — 2 and t± + t 2 = t/2 ± const, therefore the summation is 
over t 3 (or t 4 ). The summation is limited by the constraint i 3 > r on one side and by 
£3 < t — ti — t 2 « t/2 on the other. 

|£sl = 4 E E 

r<t 3 <4 <x,an,P,Pi 



P*l+1 pt2 + l p*3 + l p*4 + l 

- r (/3,/3i)^(ai,a)- r (Qi,a)^(/3i, / 3)- r (/3i,/3)^(ai,a)- r (ai,a)^(/3, / 9i) 



B 3 E E 



— y r,i = ±(t 

B 3 ^ , B 2 ^ 



"<t 3 <4 a ' ai 



(27) 
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E 2 (-) h + t 2 + 2 > t/2 
h + t 2 < t/2 
t 2 = 
t 3 + l>r 



F 2 (+) 



or 



E 3 (+) ti+t 2 + 2> t/2 
h + t 2 < t/2 
tt + Kr 
h + 2>r or 
t 3 + l>r 



ti + 2 > r 
ti < r 
t 3 < r 
t 4 = 



t 4 = 
t x + 1 > r 



t 3 + 1 < r 
t 3 + 2>r 
tx + 1 > r 



or i 2 = 





ti + 2 > r 
ii < r 
t 3 + 2>r 
t 3 <r 




Figure 12. Sets £2, £3, -F2, -D4 and F4 giving non-zero contributions. The diagram 
of Ei was "symmetrised" to highlight it's similarity with F2 



where to perform the summation over all possible bonds ((3, we invoke the probability 
conservation. The summation over all possible bonds (a, a±) is just the number of the 
bonds. 

For E 2 the situation is very similar, but we choose to sum over t 4 , 

I J? I _ A pti+l p*3+l tdU+1 pi 

1-^21 — * r (p,u)^(u 1) a) r {a,f3)^(a 1) a) r (u 1 ,a)^(P,a) r (a 1 ,a)^(afi) 

t4<|-r a,ai,P 
A 

_ \ " \ " pU+1 pi 

t4<5~r "^i./ 3 
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= ~~n2 ^(a^a)-»(a,ai)> ( 2 ^) 

where, to get to the last line, we used the symmetry of the time-reversal 
invariance, P( ai ,a)^(a,p) — P(B, a )-*(a,ai)i an d the Markov property of the probabilities. 
Unfortunately, the contribution of E 2 depends on the particular structure of short orbits 
of the graph. Hopefully, the other contributions will help: 

I p I a p*i+l r>t 2 +l pUi+l p l 

\ r 2\ — * r {a,a 1 )^{cx 1 ,a) r {a 1 , a )^{ a ,P) r { a ,l3)^{l3, a ) r {P,a)^{ a ,a 1 ) 
t 3 <r a,a lt /3 

= "r2 X] P^3)^{fi,a) P tp,a)^{a,oii) — "^2 X] ^(a\/3)->(/3,a) > (^) 

JD t 3 <ra,ai,P t 3 <r a,j3 

where to get the final result we summed over a\ and invoked the probability 
conservation. We notice that the sums in ()29|) and (J28|) are identical, up to a variable 
change, apart from the upper limit of the sum. Assuming, without loss of generality, 
that r > | — r, we obtain 

\H-m = ^ X £^SU, Q )> ( 3 °) 

f-r<t 3 <r a,/3 

but now the sum is over large values of t^, thus 

1*1-1*1 = ^ E El=^(r-Q-r)). (3D 



-r<t 3 <r a,P 



F 2 \ - \E 2 \ + \E 3 \ = — (32) 



Bringing it together with l-E^I, we obtain 

4r 
B 

The next set to evaluate is -D4, Fig. All arcs are long in this case, leading to 
the contribution 

\d \ = 2 y y p ti+i p*2+i p*3+i 

C. < / 4 (ck2,c«i)— >(o:i,Q!2) (01,02)- + (<*ii<*2) {011,0.2)— ►(c*2)£*l) 
2 



= Ji E E ^KH' (33) 

In the final contribution, coming from the set F4 of Fig. ^1 the summation is over 
£2 + ^4 ~ t — 2r. The arcs 1 and 3 are both long and of fixed length r — 1 or r — 2. This 
freedom of choice gives rise to the 4 factor in front of the sum, 

1^1 — Jr2~ X X P lltX)^{PA) P {M)^(», a iV ( 34 ) 

where the terms corresponding to the arcs 1 and 3, have already being approximated 
by 1/B. At least one of the arcs, 2 or 4 is long. Assuming, without loss of generality, 
that 2 is long, we have for the inner sum, 

Ept-z+l pU+1 _ - sr^ pU+1 _ 1 1 _ I (or\ 

a,ai,P,Pi a, ai ,P,Pi ^ 0,0! 
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Thus the result for the set F 4 is 



, , 4 ^ A(t — r — r) . . 

\ F *\ = & E 1= B 2 ■ ( 36 ) 

t 2 +U=t-2r 



The overall result is 



t 2 



2(^TR3a(r) + ^TR3b(T) + ^tr 3 c(t)) = - (\F 4 \ + \F 2 \ - \E 2 \ + \E 3 \ - 2\D 4 \) 
t 2 ( 4 , 4 4 /t \\ 2t 3 „ 

Taking into account that i^NTR3a(T) + -^NTR3b ("7") = [TU], we obtain the final result 
K TR3 (r) = 2r 3 , QED. 



5. Conclusions 

The results of this manuscript complement those of to form the derivation of the 
form factor of generic quantum graphs to the third order. 

While the derivation presented above is quite technical, the underlying idea is 
beautiful in its simplicity: the set of all orbits can be partitioned in such a way that the 
contribution of each partition set can be easily shown to be zero. And this can be done 
without evaluating the contributions of individual orbits! The nonzero "correct" result 
then arises from the intersections of the partition sets. The orbits in such intersections 
have generic properties: most of their arcs are long, making direct evaluation possible. 

We presented one of the possible partitions and obtained the expected RMT result. 
The partition is based upon three basic diagrams, each being separately partitioned. It 
is very possible that there exists a "unified" diagram with a simple partition. If such 
diagram is found, it might make possible the derivation to all orders. Alas, it has so far 
evaded all attempts to find it. 

The derivation is done for sequences of graphs satisfying two conditions: (a) graphs 
must be "uniformly ergodic" , as expressed by Eq. (JHJ) and (b) the contribution of long 
self- retracing orbits must be negligible, Eq. (fTU|) . These conditions exclude such cases 
as Neumann star graphs, which are known to possess statistics different from RMT 
predictions. The families of graphs numerically found to satisfy the RMT hypothesis, 
such as Neumann complete graphs and Fourier star graphs, also satisfy the above 
conditions. 
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(a) 




h < t/2 

h + 2a> t/2 



(b) 




a = 1 
ti < t/2 
h + 2> t/2 



Figure Al. In the first set, the tangential intersection must be of length a > 0. In 
the second set, the length of the intersection is not explicitly specified. It turns out 
that both sets contain exactly the same orbits. 




Q Q::::::Q 
oto Oy a a _! a 




Figure A2. Detailed view of the tangential intersection of length a, together with the 
adjacent vertices Sj and 
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Appendix A. Setting a = 1 in TR2 diagram 

When we evaluated the contribution of the TR2 orbits, we dropped the constraints and 
set the length of tangential intersection a to 1 (Fig. lAlj) . Here we discuss various aspects 
of this approximation. 

For a given value of a, there are 2a choices of the length t\ satisfying the inequalities 
on part (a) of Fig. IA11 For small a, the length of both arcs (ao, si) — > {fi,a ) and 
(a a , S2) — > {f'2, ota) ar e long and their contributions can be approximated by l/B. Thus 
the contribution of the TR2 diagram is 

|^E2aT a , (A.l) 

n a=l 

where T a is the contribution of the structure depicted on Fig. IA21 the tangential 
intersection of the length a. This contribution can be written as 

0-2 

T a = \ a si,<*l\ \ a fl,<xi\ II l Cr «fe,"fc+2l \ a s 2 ,a a - 2 \ Wf 2 ,a a ^ 2 \ (1 ~ S Sufl ) (1 - S S2 J 2 ) , (A. 2) 

k=0 

where the summation is taken over all indices Sj, fi and a^. Our aim is to show that 
the sum of 2aT a is equal to the corresponding contribution from part (b) of Fig. IA11 



Form factor of quantum graphs 



25 



given by 2 V] 1- Opening up the brackets in ()A.2J) and performing summation, where 



«o,"l 



i°"si,qh 


a-2 

i 4 n 




k=0 


a-2 




j 4 ni^ 


k, a k + 2 


fc=0 




1 + -Ra+2, 





possible, over Sj and /j, we get 

a— 2 a— 2 a— 2 

-^a = E IT l (T afc I afc+2 I — ^l^si.ail IT I a a k ,a k+2 I ~~ E IT I Cr Ofe>"fe+2 I I^.Qa-al (A-3) 
fc=0 k=0 k=0 

S2,a a -2\ (A-4) 

(A.5) 
where R a is 

^=EnV^ fc+ 2i 4 - (a.6) 

k=0 

The sum in (jA.lj) now folds to 

a' 

E 2aT a = 2i?o + 2a' J R a , +2 - (2a' + 2)iW (A.7) 

a=l 

Since i?o is nothing else than ^ 1, to finish the proof we need to ensure that (i) the 

remainder of (|A.7|) is negligible and (ii) we can ignore orbits with a > a'. To this end 
we take a' to be of order t and refer to condition flTUj) which is obviously sufficient. 
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